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Abstract 

During the past years the understanding of the multi-scale interaction 
problems have increased significantly. However, at present there ex- 
ists a range of different analytical models for investigating multi-scale 
interactions and hardly any specific comparisons have been performed 
among these models. In this work, two different models for the gener- 
ation of zonal flows from ion-temperature-gradient (ITG) background 
turbulence are discussed and compared. The methods used is the 
coherent mode coupling model and the wave kinetic equation model 
(WKE). It is shown that the two models give qualitatively the same 
results even though the assumption on the spectral difference is used 
in the (WKE) approach. 
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I Introduction 



The understanding of multi-scale problems is believed to be a key to gain 
deeper knowledge in basic plasma physics. The multi-scale analysis have 
had great success during the last couple of years and been able to give great 
insight in the complex plasma physical processes in e.g. the zonal flow - 
drift wave system [Ij. There are a large number of theoretical-analytical 
studies devoted to multi-scale problems using different background physics 
Refs [2]- [22] and there have also been intensive numerical investigations using 
gyrokinetic [23]- [26] and advanced fluid models [27|- [32] of the the zonal 
flow - drift wave system. 

The problem of anomalously high ion heat transport believed to be gener- 
ated by the ion-temperature-gradient (ITG) mode turbulence and controlled 
by the secondary generated zonal flows is one of the typical multi-scale prob- 
lems [2]- [1]. The zonal flows are azimuthally symmetric shear flows with 
temporal variation with time scales longer than the drift waves and having 
disparate spatial scale lengths. This is often utilized to reduce the fully cou- 
pled system to an analytically tractable problem. It is most plausible that 
robust quasi- stationary zonal flows may suppress the ion heat transport and 
give access to improved confinement regimes in tokamak plasmas. Thus, 
the understanding of the dynamical system of zonal flow - drift waves is of 
great importance for controlling the anomalous transport and the formation 
of internal transport barriers (ITB). 

The zonal flow - drift wave system has been studied previously using 
multi-scale analysis for different physics such as the electromagnetic ITG 
mode the trapped-electron-mode (TEM) [19], the electron-temperature- 
gradient (ETG) mode [TU], [13 etc. In the case of TEM background tur- 
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bulence the generation of zonal flows seem to be much weaker than for the 
electrostatic ITG mode |T9] and also for the ETG mode the zonal flow in- 
stability is quite weak. It has been speculated that the importance of zonal 
flows is of less interest in ETG turbulence. 

The problem of multi-scale interaction between a tearing mode and drift 
waves has just recently been discussed |22j . 

There are several analytical models for treating multi-scale interactions 
and among the widely used models are the coherent mode coupling method [Hj, 0- 
[TT] . the wave kinetic equation approach [5], [7]- [8], [13], [16]- [19], [22] and 
the reductive perturbation expansion method [12], [20]- [21]. There is indeed 
little comparative knowledge of the differences in the current methods used. 

The main purpose of the paper is to compare and point out differences in 
two of the main approaches to treat multi-scale problems. The comparison 
is done by deriving a model for a zonal flow generated from background ITG 
turbulence using the method of coherent mode coupling and compare this to 
an existing model derived using the Wave Kinetic Equation (WKE) p], [18] 
using the same background physics and assumptions. Direct Numerical Sim- 
ulations has been performed, showing that the scale separation assumed in 
the WKE model breaks down during the simulations [33] (and References 
therein). The method proposed here and also briefly discussed in Ref. [11] 
are complementing this previous line of research to find deeper understanding 
for the approximations made in multi-scale modeling. 

In this paper, an algebraic dispersion relation for zonal flow growth rate 
and real frequency is derived using the coherent mode coupling method in 
the presence of toroidal ITG turbulence. The model for the background 
turbulence is based on the ion-continuity and the ion-temperature-equation. 
The evolution of the zonal flow is described by a Hasegawa - Mima like 
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equation. The models employed are electrostatic and effects of trapping are 
neglected. 

It is found that the two different models for estimating the effects of multi- 
scale interaction, in this particular case the interaction of ITG driven drift 
waves and a zonal flow, are in good quantitative and qualitative agreement 
in the region of maximum linear drift wave growth rate (close to k^p = kyp ^ 
0.3). In all cases good qualitative agreement was found. 

The paper is organized as follows. In Section II the analytical model for 
the zonal flows generated from toroidal ITG modes is reviewed. Section III is 
dedicated to the results and a discussion thereof. Finally there is a summary 
in section IV. 

II Analytical models for zonal flow genera- 
tion 

The model for the toroidal ITG driven modes consists of the ion continu- 
ity and ion temperature equations p4J. In the present work the effects of 
magnetic shear, trapped particles and finite beta is neglected. It has been 
found earlier that the effect of parallel ion motion on the ITG mode is rather 
weak [51] , however, it is recently found that magnetic shear may modify the 
zonal flow generation significantly [12] . The effects of magnetic shear is out of 
the scope of the present paper. In this section the model using the coherent 
mode coupling is presented. The method how to construct the interaction of 
a zonal flow and the ITG mode turbulence and the derivation of the disper- 
sion relation for zonal flows follows earlier papers [6], [9]- [10] (and References 
therein). 
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- [0, Ti] + ^ [0, n] . (2) 

The system is closed by using quasi-neutrality with Boltzmann distributed 
electrons. The slowly varying zonal flow mode is generated from a Hasegawa- 
Mima type of equation 

- ^ V^0^^ - /iV^0zF = [0, Vl0] . (3) 

Here [A, B] = dA/dxdB/dy — dA/dydB /dx is the Poisson bracket. With 
the additional definitions n = 6n/nQ, = e6(j)/Te, Ti = 6Ti/TiQ as the nor- 
malized ion particle density, the electrostatic potential and the ion temper- 
ature, respectively. In the forthcoming equations r = Tj/Tg, v^, = psCg-y/Ln, 
Ps = Cs/Vtd where Cs = ^jTe/nii, flci = eB/niiC. We also define Lf = 
— {dlnf /dr)~^ , rji = Ln/Lx^, e„ = 2L„/_R where R is the major radius and 
ttj = r (1 + rji). The perturbed variables are normalized with the additional 
definitions n = {L^/ Ps)Sn/no, = {1^/ Pse)S(f)/Te, fi = {1^/ Ps)STi/Tio as 
the normalized ion particle density, the electrostatic potential and the ion 
temperature, respectively. The perpendicular length scale and time are nor- 
mahzed to ps and Ln/cg, respectively. The geometrical quantities are cal- 
culated in the strong ballooning limit {9 = 0, g {9 = 0, k) = 1/k (Ref. [53] ) 
where g {9) is defined hy ud {9) = uj-^eng (9)). The linear solutions to Eqs 1 
and 2 is, 

= —^^^(l-(l + ^)e^g~kl(a^ + ^Ter,g)) (4) 



2 (1 + kl) V \ 3 J ^ 3 

k I 

7 = ^^^^2 \/^^n9 iVi-Vith)- (5) 



where u = Ur + i'j and 

FLR effects in the rjuh are neglected. 

Following the same procedure as in Refs [BJ, [9]- [TU] by choosing a 
monochromatic pump wave that generates a secondary convective cell, 

^Dw{x,y,t) = ^nwe'^'^^^'^'^'-'^ + cc. (7) 
^ZFix,y,t) = ezFe^(«-^-™) + c.c. (8) 
^±{x,y,t) = ^^e*«9-^'=^)^±'=«^-^±*) + c.c. (9) 

Here ^dw is the drift pump wave {ni,(j),Ti), ^zf is the secondary generated 
wave or zonal flow (0, Tj remembering that the zonal flow does not have 
a density component), ^± are the possible sidebands {hi, 4>,Ti), {kx,ky) are 
the radial and poloidal wave numbers for the background turbulence, u is 
the complex frequency of the drift wave, is the radial wave number, Q is 
the complex frequency of the zonal flow. Where the frequencies couples as 
follows uj± = Q±uj. In Figure 1 the spectral relationship between the waves 
is shown. 

This gives a problem with four coupled waves where the zonal flow seed 
may be unstable and grow exponentially. Following the method developed in 
previous work Ref. [6] the sidebands may couple to the pump wave to give 
the zonal flow component. The sideband potentials are then found using Eq. 
1 and identifying the coupling between the drift wave component and zonal 
flow component as the sideband component. This gives the equation system, 

{iVtql- Hql)(pZF = kyql[a+(p*jj^(t)+ + a.(t)DW<P-] (10) 

0+ = i(3qxkyk±_—— (11) 

\l + a 

0_ = -^f3q.kyk^.^^^ (12) 
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Figure 1: The position of the drift wave pump {(pnw), the side bands (0+, 0_) 
and the zonal flow seed {4>zf) in wave number space. The zonal flow seed is 
further assumed to have qy = 0. 

In deriving these equations the FLR effects are considered to be small. Here 



a± = ± (13) 

k^. = kl + kl-ql (14) 

(3 = i + t + tS (15) 

6 = fa-f' + i" (16) 

+ ^re^gky 

a = Lu - ky + engky/3 (17) 

It is now possible to derive a 3rd order dispersion relation for the zonal flow 
component. 

{n + zfiql){n^ -a^) = 2(3klq,k^_[q,n-2k,a]\(PDw\^ (18) 



This third order dispersion relation is solved numerically and compared to 
the basic result for zonal flows below Ref. [TSj . 

The dispersion relation for the zonal flow derived using the wave kinetic 
model is thoroughly described in Ref. |18J and hence only the final result is 
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given here. The zonal flow dispersion relation is found using the vorticity 
equation and combining this with the wave kinetic equation as an additional 
relation between the background turbulence and the zonal flow. In describing 
the large scale plasma flow dynamics it is assumed that there is a sufficient 
spectral gap between the small scale fluctuations and the large scale flow 
i.e. for the poloidal wavenumber it is assumed that qy « ky and the wave 
kinetic equation is solved in the radial direction. It is also assumed that there 
is enough difference in the temporal variation of micro turbulence and zonal 
flow {VLzf « ^Dw)- The electrostatic potential is represented as a sum of 
fluctuating and mean quantities. In the case of zonal flows {VL) generated 
from ITG mode turbulence the dispersion relation is 

(fi + i^ql) {n - q.Vg^f = -ql{l + T + T5) \(t)DW?^ (19) 

Here Vg^ = |^ is the group velocity for the drift wave and the other param- 
eters are the same as above. An alternative statistical approach, resulting 
in a modifled wave kinetic equation, is presented in Ref. and [13] which 
also contains an extensive discussion of and comparison with the approach 
used here. 

Ill Results and discussion 

The growth rates and real frequencies from the coherent mode coupling model 
(Eq. 18) and the WKE model (Eq 19) are compared for a variety of plasma 
parameters. The dispersion relations are solved numerically and all growth 
rates and real frequencies are compared to the corresponding linear ITG 
growth rate (Eq. 5) except where explicitly stated. In this study we have 
assumed that the mode coupling saturation level (|0dvk| = YIT^) is 
reached for the background turbulence. 
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Figure 2: The zonal flow growth rate and frequency (normalized to the ITG 
mode growth) as a function of e„ with rji as a parameter for the coherent 
mode coupling model model. The parameters are t — 1, — ky — — 0.3, 
H = 0, Tji = 3.0 (square), rji = 4.0 (asterisk) rji = 5.0 (diamond). 

First in Figure 2 and 3 the zonal flow growth rate (positive sign) and real 
frequency (negative sign) as a function of e„ with rji as a parameter is shown 
for the coherent mode couphng model and the WKE model, respectively. 
The parameters are t — 1, k^ — ky — — 0.3, /j, = 0, rji — 3.0 (square), 
f]i = 4.0 (asterisk) r]i = 5.0 (diamond). In the figures it is found that the 
zonal flow growth rate is not very sensitive to a change in rji on the other hand 
in both figures large zonal fiow growth rates are found except for small e„ 
(peaked density profiles). In the models large levels of zonal fiow generation 
ilzF > Iitg) is displayed for a rather large parameter regime. Moreover, 
in the variation of e„ quite good qualitative and quantitative agreement for 
the growth rates and real frequencies in both figures are displayed. The 
quantitative agreement of in both cases are dependent on q^ c.f Figure 6, 
whereas, the qualitative agreement is conserved. 
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Figure 3: The zonal flow growth rate and frequency (normalized to the ITG 
mode growth) as a function of e„ with rji as a parameter for the WKE model. 
The parameters are t — 1, — ky — — 0.3, /i — 0, rji — 3.0 (square), 
rji = 4.0 (asterisk) rji = 5.0 (diamond). 

Second, the effect of a variation of the background turbulence spectrum 
{kx,ky) on the zonal flow growth rate (normalized to Cg/Ln) is shown in 
Figure 4 and 5. The zonal flow growth rate as a function of ky with k^ as a 
parameter is displayed for the coherent mode coupling model and the WKE 
model, respectively. The parameters arc T^j = 4, r = 1, g^; = 0.3, = and 
kx — 0.2 (asterisk), kx — 0.3 (square), kx — 0.4 (plus). The comparison is 
done for relatively small values of {kx, ky) since the perpendicular wave vector 
is considered small k\ « 1 in the derivations of the dispersion relations for 
the zonal flow growth rate and real frequency. The zonal flow growth rates 
are increasing with increasing ky (strong zonal flow generation is found if 
streamer like pump is assumed) and have a rather modest dependency on 
kx- The growth rates are decreasing for increasing kx- The zonal flow growth 
rates found for the two different models in Figure 3 and 4 are in qualitative 
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Figure 4: The zonal flow growth rate in the coherent mode couphng model 
(normalized to Cgj Ln) as a function of ky with as a parameter for rji — 4, 
T — 1, Qx — 0.3, fjL — and k^ — 0.2 (asterisk), k^ — 0.3 (square), k^ — 0.4 
(plus). 

agreement and have good quantitative agreement in the {k^, ky) region where 
the linear ITG growth rate is at maximum (close to kx = ky = 0.3). For small 
kx there is a deviation in the ky scaling of the growth rate in the coherent 
mode coupling model, whereas, for larger kx the behavior is rather similar. 

Third, the dependency of the radial wave number of the zonal flow is com- 
pared for the WKE (square) and the coherent mode coupling model (asterisk) 
in Figure 6. The parameters are r)i — A, Cn — 1.0, t — 1, kx — ky — 0.3, 
/i = 0. In the case of the WKE model the scalings of the growth rate (posi- 
tive sign) and real frequency (negative sign) are linear. The coherent mode 
coupling model has a slightly more complicated behavior where a maximum 
growth rate is found close to Qx — 0.2. The two models exhibits completely 
different scaling in Qx, however, the scalings with other parameters are con- 
served and it is only a change in the levels. In deriving the WKE dispersion 
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Figure 5: The zonal flow growth rate in the WKE model (normalized to 
Cs/L„) as a function of ky with as a parameter for rji = 4, t = 1, = 0.3, 
fj, = and kx = 0.2 (asterisk), k^ = 0.3 (square), k^ = 0.4 (plus). 

relation, a state close to marginal stability was assumed and the feature of 
a non vanishing zonal flow growth rate for small Qx may be recovered if this 
assumption is discarded. However, it should be noted that for = 0.0 the 
Reynolds stress is identically zero and the driving term disappears thus the 
zonal flow growth rate is zero at this particular value. The zonal flow growth 
rate reach a maximum if the damping term is included, in the WKE model. 

IV Summary 

The main purpose of the paper is to make a comparison of two of the current 
methods for studying multi-scale problems. In this work the two methods are 
employed to study the interaction of a zonal flow with the drift wave system. 
The first method is based on the wave kinetic equation (model derived in 
Ref. [TS]) and the other method is that of coherent mode coupling (derived 
in this paper). An algebraic dispersion relation for zonal flow growth rate 
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Figure 6: The zonal flow growth rate (normahzed to the ITG mode growth) as 
a function of for the WKE model (squares) and the coherent mode couphng 
model (asterisks) parameter for r]i — A, — 1.0, t — 1, — ky — 0.3, fj, — 0. 

and real frequency is derived using the coherent mode coupling method in 
the presence of toroidal ITG turbulence. The model for the background 
turbulence is based on the ion-continuity and the ion-temperature-equation. 
The evolution of the zonal flow is described by a Hasegawa - Mima like 
equation. The models employed are electrostatic and effects of trapping are 
neglected. 

It is found that the zonal flow growth rate is not very sensitive to a change 
in r)i. In the models large levels of zonal flow generation {'Jzf > Iitg) is 
displayed for a rather large parameter regime. Moreover, in the variation of 
e„ quite good qualitative and quantitative agreement for the growth rates for 
both models are displayed. 

The zonal flow growth rates are increasing with increasing ky (strong zonal 
flow generation is found if streamer like pump is assumed) and have a rather 
modest dependency on k^- The growth rates are decreasing for increasing 
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kx- The zonal flow growth rates found for the two different models in Figure 
3 and 4 are in qualitative agreement and have good quantitative agreement 
in the {kx, ky) region where the linear ITG growth rate is at maximum. 

It is found that the two different models for estimating the effects of multi- 
scale interaction, in this particular case the interaction of ITG driven drift 
waves and a zonal flow, are in good quantitative and qualitative agreement 
in the region of maximum linear drift wave growth rate (close to kxp — kyp ~ 
0.3). In all cases good quahtative agreement was found. 

To this end, the present work indicates that the assumption of a large 
spectral difference between the zonal flow mode and the driving background 
turbulence is of rather small importance. It seems that the important condi- 
tion is that there is enough difference in time scales of the slow varying zonal 
flow and the rapidly varying ITG turbulence. 
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